Fluctuation - induced nucleation and dynamics of the kinks on dislocation. Soliton 
and oscillation regimes in 2D Frenkel-Kontorova model. 
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Numerical simulation of the dislocation motion in 2D Frenkel - Kontorova (FK) model in the 
thermostat shows an unusual dynamical behavior. It appears that "kink" regime of dislocation 
gliding takes place in a certain region of parameters of the model but, in disagreement with the 
common views about the dislocation motion under plastic deformation condition, the kinks appear 
to be similar to sine-Gordon solitons despite the discreteness of the lattice, damping and thermal 
fluctuations. At high enough stresses and temperatures the motion of the dislocation is accompanied 
by its oscillations rather than kink nucleation. 
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I. INTRODUCTION 

The concept of solitons appeared to be very fruit- 
ful in various fields of condensed matter physics such 
as magnetism, charge-density waves, electronic structure 
of polymers, etc. Soliton is defined as a nonlinear 
wave in continuum media with particle- like properties, 
i.e. conserving its form both at the propagation and 
at the collisions with other solitons. Generally speak- 
ing, the literal transfer of the results of soliton theory to 
the case of crystal lattice is impossible since the taking 
into account the discreteness of the lattice may change 
drastically the character of solutions of nonlinear equa- 
tions. In particular, the discreteness eliminates the Gold- 
stone mode leading to the pinning of solitons, create 
the new classes of localized nonlinear oscillatory solu- 
tions (breathers), induces the coupling of "soliton" and 
"phonon" modes and so on jjj], 

Apart from this, the effects of dissipation, thermal fluc- 
tuations etc. are usually important in real situations. 
There has been recently a growing interest in non-linear 
effects in the crystal lattice dynamics particularly in the 
phenomena responsible for plastic deformation (see, e.g., 
||). According to common opinions [Q, the gliding 
of dislocations in crystals with sufficiently high Peierls 
barrier (in particular, in covalent crystals and semicon- 
ductors) is by thermofluctuation creation of kink pairs 
(kink/antikink). In the course of their motion in the 
opposite sides under the action of external stress the seg- 
ment of dislocation between the kink and antikink glides 
into the neighbor Peierls barrier valley. However, while 
the existence of such kinks is confirmed by direct elec- 
tronmicroscopic examinations pj , the information on the 
character of their motion is indirect , i.e. it is obtained 
from the comparison of results of particular theoretical 
models §, §, § with the measurements of dislocation 
mobility. The common feature of all these models is that 
the kink is considered as a solitary wave in the Frenkel- 



Kontorova (FK) model. That is, it is supposed that its 
shape would persist in the motion |[J , [|o| . This assump- 
tion, strictly speaking, has not been proven since the an- 
alytical solutions of soliton type are only available for the 
ideal continuum model, and the question on their exis- 
tence in more realistic discrete case is still open. 

In the continuum limit the ID FK model is reduced to 
the well - known sine-Gordon model jllj which is com- 
pletely integrable. This means that any initial distur- 
bance of the finite amplitude in this model breaks up into 
the succession of kinks (solitons), antikinks and breathers 
(soliton /antisoliton coupled states). In this case the soli- 
ton/ antisoliton collision may occur by the mechanism of 
elastic (i.e. without change in the shape) scattering or 
be accompanied by breather formation. It is unknown to 
what extent such a behavior might be applied to the dis- 
crete case. The existing phenomenological model of kink 
motion along a dislocation [Q, jsj are based on a quite 
contrary supposition of an unavoidable kink-antikink an- 
nihilation in the collision. This statement bases on the 
hypothesis that such factors as the discreteness of crys- 
tal lattice, damping (viscosity) and thermal fluctuations, 
so important in the actual situation, would destroy the 
features associated with the complete integrability. Of 
course, this hypothesis can be either confirmed or re- 
futed only by means of particular numerical calculations 
which, however, did not give any definite answer to this 
questions up to now |jl2|| . It should be stressed also that 
applications of results concerning soliton motion in ID 
models to the kinks are doubtful even in the continuum 
limit since the kinks are "solitons on soliton" , i.e. nonlin- 
ear excitations of 2D or 3D character. We do not know 
any analytical results about such excitations in non-one- 
dimensional continuum models and therefore computer 
simulation is the only way to investigate them. 

Ref. was the first attempt to do this but the nu- 
cleation and propagation of kinks were not observed in 
that calculation. Only recently opportunities appeared of 
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computer simulations of the dislocation motion in realis- 
tic 3D case. Rather complicated picture of the dislocation 
motion in Al has been demonstrated in However, 
since even more simple 2D situation is not understood 
now, it is worthwhile to investigate systematically pro- 
cesses of nucleation, propagation and destruction of kinks 
for a model case. 

The present work is devoted to computer simulations 
of that processes in the framework of 2D FK model. We 
show that, for certain values of the model parameters, 
the dislocation gliding is really accompanied by the nucle- 
ation and motion of kink-antikink pairs but the dynamics 
of the latter can be soliton-like. At the same time, for an- 
other region of the model parameters well-defined kinks 
are not observed at all. It means that basic assumptions 
of conventional phcnomenological models of dislocation 
gliding are at least not evident and the question about 
their applicability remains. 



II. FORMULATION OF THE MODEL AND 
SIMULATION METHOD 

Consider the common two-dimensional generalization 
of FK model (2D FK). To take into account the effects 
of damping and interaction with the thermostat we solve 
numerically (as well as in Jl2]) the Langevin equations of 
the following form 



dV 

dUni 



7 u m +£ni{t) + fi 



(1) 



where u„ is the vector of n-th atom displacement, i = 
(x, y) are the Cartesian coordinates, 7 is the damping 
factor (the atom masses being equal to 1), /, is the con- 
stant external force, £ n i(t) is the Gauss random force 
(< £ni(*) >= 0) with the correlator 



< U{t)tn>At') >= 2 1 T5 nn ,5 H ,5{t - t'), 



(2) 



As is known (TJ) such a choice of random force guar- 
antees the Gibbs distribution to be established in the 
equilibrium at temperature T. The potential energy is 
chosen as 



V = 



K 



2 £ 

<n,m> 



K 



P^^(l- C os(g.uJ), 



(3) 

where < n,m > are the pairs of nearest neighbors in 
the lattice, g is three least vectors of reciprocal lattice 
(we will consider the hexagonal lattice with a parameter 
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a = 1, then gl = ^L(-l, 0), g 2>3 = i|(i,±^)). 
first term in Eq.(3) describes the interaction of atoms in 
the harmonic approximation, and the second one is the 
simplest approximation corresponding to the sine law for 
the returning force in ID FK model. Note that, unlike 
the present work, Ref. jl2| considered the square lattice 



which, in the nearest neighbor approximation, has a cer- 
tain pathology (the zero shear modulus) . In contrast with 
this, hexagonal lattice in the nearest-neighbor approxi- 
mation represents rather general isotropic 2D situation. 
Thus, Eq. (3) can be considered as the most natural 
generalization of original ID FK model for 2D case. De- 
pending on specific interatomic interactions and chemical 
bonds, the substrate potential may be changed in com- 
parison with purely cosine form in Eq. (3). Calculating 
the substrate potentials for a series of fee metals jl5| we 
found that Eq. (3) may be considered as an adequate 
enough for substances with almost purely central pair- 
wise interactions (without strong covalency or electron 
density dependence). In principle, it is not difficult to 
take into account next terms of Fourier expansion in Eq. 
(3) but in this work we restrict ourselves only by the 
consideration of the simplest case . 

Initially lattice displacements of atoms u„(0) corre- 
sponding to the screw dislocation with Burgers vector 
b=(0,l) and the axis parallel OY, were set. The dis- 
placements for t — were chosen in accordance with the 
solution of ID FK as @ 



2 /7T 
7T V2 



— tan 



exp[ — — 



= 



(4) 



where A = — \ / 

7T V P 



Eqs.(l) were then solved for / = 0, 
T = and the "actual" equilibrium distribution of dis- 
placement field in the absence of external force was de- 
termined. After that external force f y (f x = 0) and tem- 
perature T was introduced, and the evolution of displace- 
ment field was studied. The lattice was specified the form 
of 40x40 cluster of atoms with periodic boundary condi- 
tions. Generally speaking, the latter may be dangerous 
in the simulations of dislocations in realistic models with 
inverse distance decay of deformation fields and there- 
fore with strong interactions with "mirror" dislocations. 
However, the deformation in FK model decay exponen- 
tially and long-range part of the dislocation fields is not 
taken into account. 

In order to obtain a reliable information about the evo- 
lution of the system investigated, it is necessary to use 
the so-called "strong" (mean square) methods for inte- 
grating stochastic differential equations ]l6| , [JlTj. The 
simpler and more commonly used weak methods (in par- 
ticular, those used in ]l2|| ) are actually intended only for 
the calculation of averaged characteristics such as spec- 
tral density and are not appropriate for studying individ- 
ual trajectories of the system in the thermostat. We used 
the method proposed in [ ft8[ and intended for solving the 
set of stochastic differential equations (the approach used 
in |L7| is only directly applicable for the solution of sin- 
gle equation). Specific computational formulas and the 
details of the method are presented in Appendix. 
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III. THE RESULTS OF CALCULATIONS 

As is seen from Eqs. (l)-(3) the character of solu- 
tion is determined by four independent dimensionless 
parameters 7 = j/VK, P = P/K, T = T/K, and 
/ = fy/K. Dimensionless damping 7 is of order of the 
ratio of phonon damping to phonon frequencies and is 
typically of order of 10~ 2 at room temperature and above 
|l9| . Apart from this, such choose of 7 provides strong 
attenuation of phonon waves at the lengths of order of 
the crystallite size in our simulations and eliminate es- 
sentially artifacts of periodic boundary conditions. Since 
the phonon damping in the classical region is linear in T 
(see e.g. |||) we put 7^72 = T 1 /T 2 . 

The value of the parameter P determines the tran- 
sition from continuum case to the discrete one. One 
can estimate from the expression for A in Eq. (4) that 
the width of the kink is of order of interatomic distance 
at P ~ 0.1. As it will be seen below it is roughly the 
" critical" value for the change of a character of disloca- 
tion motion. Therefore we investigate in detail the cases 
P = 0.08 and P = 0.12. As for the parameter of external 
stress it has to be compared with the Peierls stress f p . 
The latter is defined as the minimal stress which leads to 
the motion of the dislocation as a whole at T — 0. The 
results of our simulations show that this critical stress 
turns out to be temperature-dependent decreasing with 
T increase. For example, for P = 0.08, f = 0.005 the 
dislocation moves as a practically straight line already at 
/ = 0.012 — 0.016, whereas Peierls stress for that param- 
eters P and T — 0.002 was 0.030 - 0.035. The accurate 
calculation of the Peierls barrier in computer simulation 
is a complicated problem pc[| . Since we are interested 
here in the investigation of the gliding we restrict our- 
selves only the case f < f p (T) . 

Fig. [l] gives an example of the distribution of displace- 
ments around the moving dislocation. Kinks may be seen 
as inhomogeneities of the distribution along the disloca- 
tion axis. However, it is much more suitable to watch 
them displaying only dislocation line positions. The lat- 
ter is defined by the condition \u y \ = 0.5b. It corresponds 
to the center of dislocation since \u y \ varies continuously 
from to 6, see Fig. |l|. 

Fig. H demonstrates a typical (at least for the values of 
temperature and damping under consideration) picture 
of dislocation gliding for P = 0.12. It shows "soliton- 
like" behavior of kinks with their thermofluctuation nu- 
cleation, propagation without change of their shapes and 
purely elastic scattering (mutual penetration) of kink and 
antikink. To our knowledge, this is the first observation 
of such regime in 2D FK model with the interaction with 
thermostat. 

Figs. HJ^ display some pictures of the gliding for 
P = 0.08. The character of the dislocation motion for this 
value of P is much more complicated than for P = 0.12 . 



Depending on the parameters of the model three differ- 
ent regimes of the dislocation gliding are possible. First, 
kink nucleation similar to Fig. || may take place. Second, 
kinks may appear after preliminary stage of " oscillation" 
motion when the dislocation glides as a wriggled flexible 
band or, by another words, kinks arise as a result of evo- 
lution of the oscillations of the dislocation segment (Fig. 
|3|) . The most important feature here is the demonstration 
of a possibility of a third regime, namely, rather fast glid- 
ing without participation of kinks at all, only due to the 
oscillations of the dislocation line (Fig. ^). Sometimes 
the average dislocation velocity without participation of 
kinks is even higher (compare the motion of dislocation 
line in Fig. ^| where kinks exist and in Fig. ^ without 
kinks). Since the oscillating dislocation does not lie in 
a single Peierls valley the "effective" Peierls barrier for 
such dislocation segment turns out to be lower than for 
straight one. Results obtained here shows that the regime 
of " wriggled flexible band" occurs at large enough tem- 
perature and external stress close to the Peierls stress. 
Similar gliding mechanism due to oscillations has been 
observed recently in molecular dynamics simulations of 
the 3D dislocation motion in Al [HJ. According to our 
simulations, the transitions between regimes occur at the 
change of temperature and damping, e.g., for T = 0.002 
three regimes mentioned above take place for 0.02 < 7 
< 0.03, 0.03 < 7 < 0.04, and 7 > 0.04, correspondingly. 
General trend is the change of kink regime by the oscil- 
lation one at the temperature and/or damping increase. 
The value of external stress is less important provided 
that it is lower than the Peierls stress for a given tem- 
perature. 

Another important feature of that regime is the exis- 
tence of long-lived coupled kink-antikink pairs similar to 
discrete breathers 0. An example of such behavior is 
seen in the part of Fig. || corresponding to 24 < t < 34. 
It is worthwhile to stress that we mean breather-like ex- 
citations on the " soliton" (dislocation) and not standard 
discrete breathers as more or less isotropic localized os- 
cillating in time solution of non-linear equations. Sim- 
ilar breather-like excitations may be observed also for 
P = 0.12 but their lifetime is much shorter than for 
P = 0.08. Note that the lifetime increases with 7 increase 
for given values of other parameters of the model. 

We never observed the kink-antikink annihilation in 
our simulations; if kinks appear they always behave as 
solitons at the collisions. It is rather unexpected in com- 
parison with common views Q], |q], [jn]. To clarify the 
situation we investigate kink-antikink collisions generat- 
ing kink-antikink pairs right in the initial configuration. 
The results of the simulations shown in Fig. ^| demon- 
strate that both penetration and annihilation are possible 
depending on the parameters of the model, the annihila- 
tion being take place at rather high values of the damp- 
ing and temperature (Fig. |^b). However, it is that the 
region of the parameters when, instead of thermofluctua- 
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tion nucleation of kinks, oscillation regime of the gliding 
takes place. Therefore we believe that the kink-antikink 
annihilation is not important process for the dislocation 
gliding under "normal" conditions when the creation of 
kinks is of purely thermofluctuation nature. However, 
it may be important if kinks could create, e.g., by colli- 
sions of the dislocation with another defects, etc. This 
question needs further investigations 



purities M might be also studied be means of similar 
simulation. At the same time, we believe that the results 
obtained may be of more general interest than for physics 
of plasticity only since until now rather incomplete infor- 
mation about soliton-likc excitations in two-dimensional 
discrete systems are available. 
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IV. DISCUSSION AND CONCLUSIONS 

Summing up the results obtained, it might be stated 
that the physical picture of the dislocation gliding ap- 
peared to be much more complicated that it is supposed 
in traditional models of dislocation motion in the Peierls 
relief (§, §]. First of all, we observe (for P = 0.12) ther- 
mofluctuation nucleation of kinks supposed in that mod- 
els but, in contrast with them, the character of motion of 
the kinks appears to be purely " solitonic" . It is nontriv- 
ial since for these values of the parameters the widths of 
kinks are of order of the lattice constant and the results 
of continuum models such as sine-Gordon model are def- 
initely not applicable. This implies the annihilation of 
kinks and antikinks is unlikely in this regime, and their 
lifetime is anomalously long. It means that under these 
conditions the thermodynamical equilibrium in the sys- 
tem of kinks and antikinks does not practically establish, 
and their concentration should be essentially higher than 
that predicted by the phenomenological models. Note 
that in the special pulse regime of loading the anoma- 
lously high concentration of kinks was observed in Si ]2l| ] 
and Ge S. The explanation proposed in || is based on 
the essential effect of impurities. According to the results 
obtained in the present work this explanation is not the 
only possible one and high concentration of kinks may be 
intrinsic property of that materials. 

For larger values of width of kinks (P = 0.08) the 
character of dislocation motion is very complicated. In 
contrast with the case P = 0.12 the creation of kinks in 
this case is not typical thermofluctuation nucleation but 
rather a result of more or less long evolution of oscilla- 
tions of the dislocation segment. Both kinks and oscil- 
lations of dislocation segments play an important role in 
the gliding. It is shown that the gliding with rather high 
velocity may be provided by oscillations only without 
participation of kinks ("wriggled flexible band" regime). 
More detailed investigation of these processes depending 
on the temperature, damping, external stress and the 
shape of substrate potential relief would allow to correct 
and improve existing phenomenological models of plastic 
deformation. 

In the conclusion it should be noted that the methods 
used here can be applied for studying variable loading 
conditions with a time-dependent external force (ultra- 
sound, pulse loads). The question on the role of im- 
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VI. APPENDIX 

Let us represent Eq. (1) as a set of first -order equa- 
tions in a form used commonly in the theory of stochastic 
differential equations (SDE) 



dv n 



du ni = v ni dt 
(-V m V r (u) - jv ni )dt + nduj n i(t) 



(Al) 



where K — ^/2^T, dw n i(t) is a standard Wiener process. 
Introducing Anm -dimensional vectors 



v 

-W(u) 



7V 





b. 



b ri — n5ri, r=l,...,q, 1=1, ...,q, q = 2 ram, 
one can rewrite Eqs. (Al) in a common form as 

i 

dl? = qf(f, -&)dt + ojdu r (t) (A2) 

r=l 

To integrate SDE's (A2) numerically we use the mean- 
square Runge-Kutta type method of the 3/2- order which 
leads to the following approximate solution 18 



2 i 



r=l 



r—1 r—1 



4 



, 9 r Ioor)k +y^(A r £i"a > / r oo)fc + y^(£iA r "a > /o r o)fc. 



r=l 



r=l 



Here 16 k is the mean-square approximation of the vec- 
tor ~x* at the instant 

i?i = h~ct{t k , lt k ), i? 2 = hlt(t k + h/2, lt k + k[/2), 
t 3 = hlt(t k + h/2, x> fe + lt 2 /2), 
I? 4 = h~ct(t k+ i,l? k + K 3 ) 



d 9 d 1 i p g. 



2 ^ ^ ~r rdxWxi 
r—1 — l 



lor — hl r — I r Q, I r QQ — hI r Q — J r 



lorO — %Jr — hI r Q, IoQ r — ~^^r ~ Jn 



3 4V3 12^5/ ' 

h = t k +i — t k ] p — Anm is the dimension of ~af vector; 
^ ri ri r and t; r are independent normally distributed ran- 
dom variables with zero mean value and dispersion equal 
to unity. 

Generally speaking, this method is "ft 3 / 2 — method", 
i.e. it provides the accuracy of order of /i 3 / 2 for H? ' k for 
every realization of stochastic process. However it may 
be proven |l8| that for the case of small enough noise it 
is even more accurate (" ft 4 method" ) . We have chosen in 
our simulations h = 0.005; in that case the accuracy of 
the calculations for trajectories as well as of Gibbs distri- 
bution at the equilibrium appeared to be restricted only 
by Monte-Carlo inaccuracy. 
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CAPTIONS TO FIGURES 



FIG. 1. A typical picture of 

the displacement distribution for the moving dislocation at 
P = 0.08, T = 2 ■ 10 _3 ,7 = 0.02,/ = -0.008. Coordinates 
x,y are in units of a, u y (vertical axis) is in units of Burgers 
vector. 



[1] Soliton in Action, Ed. by K. Lonngren and A. Scott 
(Academic Press, New York, 1978); K. Maki, In: Elec- 
tronic Properties of Inorganic Quasi- one- dimensional 
Compounds, Ed. by P. Monseau (Academic Publishers, 
Kluwer, Dordrecht, 1985), p. 125. 



FIG. 2. The configuration of the dislocation line is shown 
at different times t = 10, 13, 14, 15, 16, 18, 20, 22^ 25, 27, 29, 
31, 33, 35, for P = 0.12,T = 2 ■ 10" 3 ,7 = 0.02,/ = -0.151. 
The pictures correspond to the consequent events from left 
to right, the second row is the continuation of the first one. 
Dashed line shows the initial position of the dislocation. Cir- 
cles indicate position of the atoms. 
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FIG. 3. As in Fig. 2, t= 10, 13^16, 18, 20, 22, 24, 26, 
28, 30, 32, 34, 37, 40, for P = 0.08, T = 2 ■ 10~ 3 , 7 = 0.035, 
/ = -0.008. 



FIG. 4. As in Fig. 2, t = 7, 10, 13, 16, 19, 22, 25, 28, 31, 34, 
37, 40, 43, 46, for P = 0.08, T = 2-10~ 3 , 7 = 0.045, / = -0.01. 

FIG. 5. The dynamics of the kink-antikink pair created in 
the initial configurations. The dislocation line at different 
times t =,4.5, 6, 7.5, 9, 10.5, 12, 13.5^15 is shown. Here 
P = 0.08, T = 2 ■ 10~ 3 and (a) 7 = 0.02, / = -0.030 (penetra- 
tion regime) (b) 7 = 0.03, / = —0.040 (annihilation regime). 
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